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The tunneling conductance in the half-metal/conical magnet/superconductor (HM/CM/SC) is 
investigated by the use of the combined Blonder-Tinkham-Klapwijk (BTK) formalism and the 
Bogoliubov-de Gennes (BdG) equations. We show that the conductance calculated self-consistently 
differs significantly from the one calculated in the non-self-consistent framework. The use of the 
self-consistent procedure ensures that the charge conservation is satisfied. Due to the spin band 
separation in the HM, the conductance in the subgap region is mainly determined by the anomalous 
Andreev reflection the probability of which strongly depends on the spin transmission in the GM 
layer. We show that the spin of electron injected from the HM can be transmitted through the GM 
to the SG adiabatically or non-adiabatically depending on the period of the exchange field modu¬ 
lation. We find that the conductance in the subgap region oscillates as a function of the GM layer 
thickness wherein the oscillations transform from irregular, in the non-adiabatic regime, to regular 
in the adiabatic case. In the non-adiabatic regime the decrease of the exchange held amplitude in the 
GM leads to the emergence of the conductance peak for one particular GM thickness in agreement 
with experiment [J.W.A Robinson, J. D. S Witt and M. G. Blamire, Science 329, 5987]. For both 
transport regimes the conductance is analyzed over a broad range of parameters determining the 
spiral magnetization in the GM. 


I. INTRODUCTION 

In recent years the quantum transport in the ferromag- 
net/superconductor (FM/SC) junctions has attracted 
growing interest due to a possible existence of the spin- 
triplet pairingir— and novel transport phenomena re¬ 
lated to this unique superconducting state^^i^ For the 
s-wave superconductors with the spatially symmetric or¬ 
bital part of the Cooper pair wave function, the Pauli 
principle requires that its spin part is antisymmetric, 
which means that the spin-singlet seems to be the only 
possible state for the Cooper pair. However, many years 
ago Berezinskii proposed a possible existence of the spin- 
triplet state in a system with the s-wave interaction which 
do not violated the Pauli principle. The triplet pair¬ 
ing correlations proposed by Berezinskii are odd in time 
(or frequency) and can appear in systems with sort of 
the time-reversal symmetry breaking mechanism. Recent 
studies suggest that the spin-triplet Cooper pair correla¬ 
tions can be induced and observed experimentally in the 
FM/SC junctions with the spin-active or magnetically 
inhomogeneous interface A”— 

In the normal metal/superconductor (NM/SC) junc¬ 
tions the electrons incident on the interface from the NM 
side are reflected as holes with opposite spins. This mech¬ 
anism, known as the Andreev reflection^ leads to the 
proximity effect - the superconducting pairing correla¬ 
tions penetrate into the normal metal over the distance 
as long as one micron at low temperatureThe prox¬ 
imity effect significantly changes if we replace the normal 
metal by the ferromagnet. The exchange interaction in 
the ferromagnet results in the different Fermi wave vec¬ 
tors for electrons with opposite spins forming the Cooper 


pairs. This wave vector mismatch is compensated by 
the non-zero total momentum of the electron pairs giving 
rise to the oscillations of the spin-singlet superconducting 
correlations in the ferromagnet^!^ known as the FFLO 
oscillationsSince the exchange field tends to align 
the electronic spins along the field direction, the spin- 
singlet superconducting correlations in the ferromagnet 
are strongly suppressed leading to the short-range pene¬ 
tration length. In contrast to the short-range proximity 
effect for the spin-singlet state, the spin-triplet state with 
m = ±1, with both electronic spin of the Cooper pair di¬ 
rected along the exchange field, is robust against the pair 
breaking induced by the exchange interaction. Therefore, 
the spin-triplet superconducting correlations {m = ±1), 
if they exist, can penetrate the ferromagnet over the dis¬ 
tance comparable to this observed in the NM/SC junc¬ 
tions. This phenomenon, called long-range proximity ef¬ 
fect was predicted theoretically by Bergeret et al. (see 
Refs, d, [3, and[i). It was found^i^i^ that in the FM/SC 
multilayer junctions with the spin-active or magnetically 
inhomogeneous interface (the spin-flip processes are pos¬ 
sible) all three components m = 0 and m = ±1 of the 
spin-triplet state can arise. Despite few theoretical stud¬ 
ies on the spin-triplet pairing induced by the spin-active 
interface, including the effect of domain wall^ spin-orbit 
couplingii or spin-dependent potential^ up to date, the 
direct evidence of the spin-triplet supercurrent has been 
observed in multilayer FM/FM/SC systems with a non- 
collinear magnetization of the ferromagnetic layers 

A first experimental hint for the long-range proxim¬ 
ity effect was reported in a half-metal Josephson junc¬ 
tion based on Cr02'^ However, since the measured crit¬ 
ical current varied by two orders of magnitude in simi- 


2 


lar samples, the results of this experiment needed to be 
confirmed. The strong evidence for the long range prox¬ 
imity effect was then reported in the Josephson junc¬ 
tions based on Co^^ The dependence of the critical 
current on the Co layer thickness, which agrees with 
the theoretical expectations, provides a strong experi¬ 
mental confirmation of the existence of the spin-triplet 
pairing in the FM/SC heterojunctions. Further stud¬ 
ies on the spin-triplet pairing concerned the FM/SC/FM 
and FM/FM/SC junctions with a relative magnetiza¬ 
tion between the ferromagnets. The spin triplet pairing 
in the clean FM/SC/FM nanostructures with an arbi¬ 
trary angle between the magnetization of the FM lay¬ 
ers was theoretically studied by Haltermann et al. in 
Refs. [ 234 ^ The authors used the self-consistent solu¬ 
tions of the microscopic Bogoliubov de-Gennes (BdG) 
equations and analyzed the spin-triplet correlations as a 
function of the relative magnetization between the mag¬ 
nets. The self-consistent calculations allowed to con¬ 
firm the experimentally observed angular dependence 
of the critical temperature Tc which monotonically in¬ 
creases due to the presence of the long range spin- 
triplet correlations. Tc reaches minimum if the relative 
magnetization is parallel and maximum for antiparal¬ 
lel magnetization^2i A different behavior was observed 
for the FM/FM/SG nanostructures for which the crit¬ 
ical temperature is minimized in case of perpendicular 
alignment of the magnet iz at ion 

Research on the spin-triplet pairing in the FM/FM/SG 
junctions has been recently extended to systems with 
the conical (helical) ferromagnets (GM). Efforts to 
control the long-range triplet supercurrent has been 
recently demonstrated in Josephson junctions based 
on holmium(Ho)-Gobalt(Co)-holmium(Ho) multilayer 
setup 4^ One has been observed a nonmonotonic depen¬ 
dence of the critical supercurrent as a function of the 
Ho layer thickness, dcM, with peaks for dcM =4.5 nm 
and 10 nm. By increasing the Go layer thickness a 
slow decay of the critical current has been reported in 
agreement with theoretical calculationsNevertheless, 
the theoretical model presented in Ref. [3l| does not ex¬ 
plain the complex dependence of the critical supercur¬ 
rent on the Ho thickness. The nonmonotonic behavior of 
Ic{dcM) has been obtained by Halasz et al. in Ref. [32| 
who have performed calculations in the clean limit us¬ 
ing Eliashberg equations. Similar dependence has been 
also demonstrated by the use of the Blouder-Tinkham- 
Klapwijk (BTK) approach^^i^ 

In the mentioned theoretical works^~— the proximity 
effect at the GM/SG interface has been neglected mean¬ 
ing that the superconducting pair potential has been as¬ 
sumed to be a step function. However, as shown by recent 
studies^ only the self-consistent calculations of the tun¬ 
neling conductance guarantees that the charge conserva¬ 
tion law is satisfied. It means that one cannot properly 
determine the tunneling conductance in the EM/GM/SG 
heterostructures by using the non-self-consistent frame¬ 
work. The full self-consistent approach is needed. The 


self-consistent calculations of the spin-triplet correlations 
in two layered CM/SC junctions have been presented in 
Refs. [33 and [ 33 . Nevertheless, these studies concern only 
the spin-triplet correlations between the CM and SC. 
They do not include the analysis of the tunneling con¬ 
ductance (transport calculations), the influence of the 
EM layer attached to the CM or the influence of the CM 
layer thickness. Summing up, the theoretical analysis of 
the tunneling conductance through the EM/CM/SC het¬ 
ero junctions with the inclusion of the proximity effect in 
the full self-consistent framework has not been presented 
until now. 

In the present paper we report the full self- 
consistent calculations of the tunneling conductance in 
the HM/CM/SC junctions. The charge transport in the 
considered system is mainly determined by the anoma¬ 
lous Andreev reflection, the probability of which strongly 
depends on the spin transmission in the CM layer. We 
consider the conductance in two cases in which the spin 
transport is adiabatic and non-adiabatic. The conduc¬ 
tance is analyzed over a broad range of parameters de¬ 
termining the spiral magnetization in the CM. We show 
that the tunneling conductance in the HM/CM/SC junc¬ 
tions strongly depends on the spin transport regime. The 
paper is organized as follows: in Sec. [IT] we introduce the 
basic concepts of the theoretical scheme based on the self- 
consistent solution of the BdG equations and the BTK 
formalism. In Sec. ini we present the results while the 
summary is included in Sec. IIVI 


II. THEORETICAL METHOD 

We consider the EM/CM/SC structure schematically 
illustrated in Eig. [H In the x — z plane the system is 
assumed to be infinite while the y axis is perpendicular 
to the layers whose lengths are denoted by dpu^ dcM, 
dsci respectively. The value and the direction of the 



FIG. 1. (a) (Color online) Schematic of the FM/CM/SC 

nanostructure. The exchange field h, denoted by the red ar¬ 
row, is directed along the 2 : axis in the FM and has a spiral 
structure in the CM. (b) The conical magnetic configuration 
with the cone angle a and a rotational angle /3. a is the lattice 
constant. 

exchange field h, denoted by the red arrow in Fig. [U 
depends on the position. It is directed along the z axis 
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in the ferromagnet, h = ( 0 , 0 , while, in the conical 

ferromagnet, h is given by 

{ hcM sin a sin (/3o + 

heM cos a , ( 1 ) 

hcM sin a cos (/3o + 

where hcM is the exchange field amplitude, /3o deter¬ 
mines the angle of a relative magnetization of the CM 
layer measured from the FM layer at the FM/CM in¬ 
terface while a and /3 are the cone and rotation angle 
whose physical meaning is depicted in Fig. [IJb). From 
Eq. m the spatial period of the helix exchange field is 
A = 27ralP, where a is the lattice constant. 

In the present paper we consider the charge transport 
through the FM/CM/SC junction and analyze the tun¬ 
neling conductance as a function of the system parame¬ 
ters. As mentioned above, the correct analysis of the con¬ 
ductance behavior require the inclusion of the proximity 
effect. This can be done only by the full self-consistent 
calculations in which the pair potential distribution is de¬ 
termined from the microscopic BdG equations. The self- 
consistent procedure used in the paper consists of two 
steps. First, we determine the self-consistent pair poten¬ 
tial A{y) in the nanostructure by solving BdG equations. 
Then, A{y) is used to calculate the tunneling conduc¬ 
tance within the BTK approach^^"— Below, both these 
steps are described in detail. 

A. Self-consistent pair potential A{y) calculations 


where Tj(r), Ts(r) are the creation and annihilation op¬ 
erators with spin 5 , h(r) = (ha,(r), h^(r), hz{T)^ is the ex¬ 
change field, (T = {(Tx^cFy^ (Jz) is the Pauli matrices vector, 
/i(r) is the chemical potential and A(r) is the spin-singlet 
pair potential in real space defined as 


A(r) = 5 ((r) , (3) 


where ^(r) is the phonon-mediated electron-electron cou¬ 
pling constant. The generalized Bogoliubov transforma¬ 
tion 


^^(r) = [Uns{r)ln + VsKsi^hl] > ( 4 ) 

n 

where 7 ^ and 7 J are the quasiparticle annihilation and 
creation operators, r4ns(i*) and are the electron and 

hole components of the amplitudes vector and r]s = 1 (~ 1 ) 
for spin down(up), reduces the Hamiltonian (|2]) into the 
diagonal form. By the commutation relation 


[T^, n] = {He - hz)^t - {hx - ihy)% + AT;, (5) 


The effective BGS Hamiltonian of the considered sys¬ 
tem is given by 

'H = Yj - n{r)j ^s(r) 


- f dVy]4'^(r)[h(r) 

ss' 



[4'4, n] = {He + - {h^ + ihy)^i, + A4't, (6) 


where Hg = ——/Li(r) is the single electron Hamilto- 
nian, and using the fact that the system is infinite in the 
X — z plane we obtain the BdG equations in the quasi-one 
dimensional form 



He-hM 

-hx{y) + ihy{y) 

0 
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0 
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0 

A*(y) 
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A*(y) 

0 

-hx{y) + ihy{y) 

-He-h,{y) ) 
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Untiv) \ 

Uniiv) 

Vntiy) 

\ ^nliy) / 


(7) 


Equations 0 are coupled with the expression for the pair 
potential given by 


A(y) 


aiv) 


E 

\En I <huJD 


Unt{y)vli{y) + Wn4.(yXt(2/)] 


X [1 - 2/(E„)], ( 8 ) 


where f{E) is the Fermi-Dirac distribution. The 
summation in Eq. (jS]) is carried out only over the elec¬ 
tronic states with energy En inside the Debye window 
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\Ej^\ < huoDi where ujd the Debye frequency. In our 
approach ^(r) is assumed to be nonzero only in the SC 
layer. 


The self-consistent procedure used to solve the BdG 
equations © is similar to these reported in previous 
papersThe main difference is that the assumed 
basis functions have the form of the plane waves. Such 
choiceis needed in the transport calculations because it 
guarantees nonzero current at the boundaries of the sys¬ 
tem. The self-consistent procedure can be described as 
follows. First, the BdG equations © are diagonalized in 
the basis of the plane waves 


Untiy) 
Unliy) 
Vn-tiy) 
\ V„l{y) 



'^nqX 


exp{ikqy) 


(9) 


where Unq^^ Unqi, Vnq^^Vnqi are the expansion coefficients, 
kq = 2'KqjL is the wave vector with q being an integer 
while L is the total length of the nanostructure. Then, 
using the calculated wave functions 
the new pair potential A(^) is determined on the basis of 
Eq. dD). This new A(^) distribution is used in the next 
iteration in which we again solve the BdG equations and 
determine A(^). This procedure is repeated until the 
convergence is reached. Due to the high computational 
complexity of such scheme the parallel implementation 
of the numerical procedure is required. 

Finally, we also calculate the magnetization vector m 
given by the formula 

mx{y) = < {un^{y)u*^^{y)+Uni{y)ul,^{y))f{En) 

n ^ 

- “ /(-£^n)]|aO) 


'my{y) = -iys \ (Mnt(yXi(2/) - Uni{y)K^{y))f{En) 

n ^ 

+ Kt(y)X(2/) - Vni{y)vl^{y))[l - /(-B„)]|, (11) 


rriziy) = ^bY] Xt(yXt(y) “ Uni{y)u*„^{y))f{En) 

n ^ 

+ (vnt(yXt(y) - 'Vni{y)v*n^{y))[i - /(X)] j()-2) 


where are the Bohr magneton. 


B. Tunneling conductance calculations 


on the grid yy = va with lattice constant a {v = 
1,2,...). We introduce the discrete representation of 
the quasi-particle wave vector as follows: = 

{W{yu)), \u^{yu)),\v^{yu)), \v^{yu)))^ = Introduc- 

ing a set p of Pauli-like matrices in electron-hole space, 
the discretized tight-binding form of the Hamiltonian in 
Eq.([7j) is given by 

H = Y,\pz® [(2^ - - hzuUz] 

V ^ ' 

-^{p.0 (tl)|^'.+i)(^.|+F.c} (13) 

V ^ ^ 

T ^ ^ I [ i^^l/Px l) G) {hyvPz) G CFy 

where py = p{yv), {hxu.hyy.h^u) = 
[hx{yv),hy{yy),hz{yy)), t = and 1 is 

the unity matrix. 

Let us assume that the electron with spin-up is in¬ 
jected from the EM into the SG through the CM layer. 
There are five possible scattering processes: normal re¬ 
flection with spin conservation normal reflection 

with spin-flip (RH), reflection as a hole with opposite 
spin (normal Andreev reflection, reflection as a hole 

with spin conservation (anomalous Andreev reflection, 
RJJq) and transmission as a quasi-particle In the 

above, denotes the reflection probability where 

upper and lower right index corresponds to the state of 
an incident particle while upper and lower left index is 
associated with the reflected one, is the transmis¬ 
sion probability where upper indexes indicate the state 
of incident particle. Analogous scattering processes can 
be distinguished for the spin-down electron injected from 
the EM. Their probabilities are marked by 
and , respectively. 

According to the BTK approach the current through the 
FM/GM/SG junction can be calculated from the formula 

liv) = I G{E) [f{E - eV) - f{E)] dE, (14) 


where V is the bias voltage and f{E) is the Fermi-Dirac 
distribution. At low temperature the energy dependent 
tunneling conductance G{E) = dl/dV\v=E (in units of 
e^/h) is given by 


G{E) = 


1 + hpM 


1 — hpM 


l+RlliE) + Rj^^iE) 
-RtliE) - Rj^^iE) 
1 + RiiiE) + RliiE) 


The tunneling conductance calculations have been per¬ 
formed within the tight-binding approximation using 
the Kwant packageFor this purpose we have trans¬ 
formed the BdG equations © into the discretized form 


-RttiE) - Rli(E), 


(15) 


where hpM^ expressed in units of //, corresponds to the 
spin polarization at the Fermi level in the FM layer. In 
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our calculations the reflection probabilities in Eq. (USD 
are determined by the use of the Kwant package^ which 
requires the implementation of the discretized tight- 
binding Hamiltonian given by Eq. ([T3|) . In the paper, 
we consider the forward tunneling conductance with the 
angle of the incident electron = 0 and neglect the scat¬ 
tering potential at the interfaces. 


III. RESULTS AND DISCUSSION 

In this section we analyze the tunneling conductance 
through the EM/CM/SC junctions by the use of the full 
self-consistent approach presented in Sec. [TTl Since the 
first experimental evidence for the long-range proximity 
effect (spin-triplet pairing) was reported in a half-metal 
Josephson junction^ we restrict our analysis to the case 
in which the ferromagnetic layer is embedded in a half¬ 
metal (HM), hpM = 1. In our calculations we neglect the 
Eermi wave-vector mismatch between the layers assum¬ 
ing a constant value of the chemical potential fi through¬ 
out the nanostructure. Its value is used as the energy 
unit. In the calculations we adopt the following values 
of the parameters: zero temperature energy gap in the 
bulk Ao = 0.01, Debye energy Tiujjj = 0.1, temperature 
ki)T « 10“^ and the lattice constant a = 0.35 nm cor¬ 
responding to the conical magnet Holmium^S Other pa¬ 
rameters determining the magnetic configuration of the 
MC, namely the exchange field amplitude hcM^ the cone 
angle a, the rotation angles (3 and /3o, as well as the 
CM layer thickness dcM are used to analyze the tunnel¬ 
ing conductance through the HM/CM/SC junctions and 
vary from one simulation to another. 

As predicted by Bergeret et the join effects of the 
Andreev reflection and the proximity at the EM/SC in¬ 
terface allow for the coexistence of the spin-singlet pair¬ 
ing correlations (| ti) “ I ti))/\/2 and the spin-triplet 
pairing correlations with the total spin projection m = 0 
(I ti) + I ti))/v^)- If a magnetically inhomogeneous 
layer, such as the CM layer, is present between the EM 
and SC, the spin-triplet state with m = 0 can be rotated 
to the state with m = 1 (| tt))- H means that the ex¬ 
istence of the spin-triplet pairing with m = 1 strongly 
depends on the spin transmission in the CM layer. This, 
in turn, is determined by the exchange field which in the 
CM has a rotating component varying with the period 
A. Depending on A the spins of electrons injected from 
the EM can be transmitted through the CM adiabatically 
- the spin orientation follows the spatial modulation of 
the exchange field, or non-adiabatically - the period of 
the exchange field modulation is so short that the elec¬ 
tron spin is not able to adopt to the field changes. The 
degree of adiabaticity can be defined by the parameter 
Q = where cj/j, = is the magnetic field 

modulation frequency in the electron’s frame of refer¬ 
ence, Vf is the Eermi velocity and ujl = hcM/ki is the 
frequency of the spin Larmor precession. In the adia¬ 
batic regime, Q >> 1. Below, we analyze the tunneling 


conductance through the HM/CM/SC junctions in both 
adiabatic and non-adiabatic regimes. 

A. Non-adiabatic regime 

All results presented in this subsection have been ob¬ 
tained for the rotational angle /3 = 30^ which corre¬ 
sponds to the spatial period of the helical exchange field 
A = 3.4 nm measured in the Holmium^^ Eor this value 
of A the spin transport through the considered system is 
non-adiabatic, Q < 1. 



y(nm) 



y(nm) 

FIG. 2. (Color online) (a) Comparison of the self-consistent 
(red solid line) and the non-self-consistent (blue dashed line) 
pair potential A(y). (b) mx{y) and (c) mz{y) components of 
the magnetization in the nanostructure. For comparison the 
exchange field components hx{y) and hz{y) are presented by 
the blue dashed lines. In panels (b) and (c) the vertical gray 
dotted lines mark the boundaries of the CM layer. 

In Fig. [2] we present the pair potential and the mag¬ 
netization in the nanostructure calculated for hcM = 1, 
a = 90° and the thickness of the conical magnet layer 
dcM = 20 nm. As shown in Fig. [2](a) the self-consistent 
pair potential significantly differs from the one used in 
the non-self-consistent approach. Due to the proximity 
effect A{y) does not have the step-like form but smoothly 
increases in the SC region reaching its bulk value Aq for 
a distance grater than the coherence length. Similarly, 
as the magnetism alters the superconductivity near the 
CM/SC interface, the superconductivity also influences 
the magnetism. This so-called reverse proximity effect 
allows to penetrate the magnetization into the SC region 
as presented in Fig.[2jb,c). Almost a tenfold reduction of 
the magnetization amplitude in the CM layer (red lines), 
as compared to the exchange field (blue dashed lines), 
results from the fact that the chosen value of A corre¬ 
sponds to the non-adiabatic transport regime. In this 
regime the changes of the exchange field seen by elec- 
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trons flowing through the nanostructure are so fast that 
the their spin do not have enough time to adopt to these 
changes. As a result, the electron spin rotates around the 
exchange field irregularly. Note that, in accordance with 
Eqs. the magnetization is expressed as a sum of 

the averaged spin over states with different wave vectors. 
Since spins of these states rotate around h with differ¬ 
ent irregular frequency, this sum averages to low value. 
In subsection B we will show that the suppression of the 
magnetization does not dependent on the conical magnet 
thickness dcM but, as expected, is mainly determined by 
the spatial period of the helical exchange field A. 



FIG. 3. (Color online) Normalized conductance G as a func¬ 
tion of energy E (in units of Ao) for different thicknesses of the 
CM layer dcM = 1,1.5 and 2 A. Results obtained by the use 
of the self-consistent (solid lines) and the non-self-consistent 
(dashed lines) approach. The conductance calculated with¬ 
out the CM layer is marked by the gray lines. Results for 
hcM — a — 90"" and /3o = 0. 

For the self-consistent pair potential A(^) we calculate 
the tunneling conductance using the procedure described 
in Sec. in B. Figure |3] presents the normalized conduc¬ 
tance G as a function of energy for different thicknesses 
of the CM layer calculated by the use of the non-self- 
consistent (dashed lines) and the self-consistent (solid 
lines) procedure. For comparison, we also mark the con¬ 
ductance calculated without the CM (gray lines). As 
we see the dashed and solid gray lines overlap which 
results from the fact that the conductance in this case 
is nonzero only for the high-energy limit, above the en¬ 
ergy gap, for which the electron incident into the super¬ 
conductor does not experience much difference between 
the step-like pair potential and the smooth pair potential 
from the self-consistent approach. Results presented in 
Fig. [3] clearly show that the self-consistent conductance is 
considerably different than this obtained in the non-self- 
consistent framework. The most pronounced difference 
between them is observed in the subgap energy range. 
Based on the results presented in Fig. [2] and [3] one can 
formulate the following conclusion: to properly deter¬ 
mine the tunneling eonduetanee in the HM/CM/SC het- 
erojunetions the full self-eonsistent ealeulations ineluding 
the proximity effeet are needed. 

To explain the conductance behavior presented in 



E/Ao E/Aq 

FIG. 4. (Color online) The reflection and transmission prob¬ 
abilities RJJ, RH, ^ function of energy (in 

units of Ao) for the CM layer thicknesses (a) dcM — A and 
(b) dcM = 1.5 A for which the conductance in the subgap re¬ 
gion reaches the minimum and maximum value, respectively 
(see Fig. [3|). Results for hcM = 1, a = QO"" and /3o = 0. 


Fig. [3] let us first discuss the half-metal/superconductor 
(HM/SC) junctions without the CM layer. In the 
HM/SC junctions the normal Andreev reflections are for¬ 
bidden due to the isolation of the spin band. This leads 
to the zero conductance in the subgap energy regime as 
depicted by the gray solid lines in Fig. [3] (we assume no 
scattering potential at the interface). The situation dia¬ 
metrically changes if we put the conical magnet between 
the HM and SC. As predicted by Bergeret et al^^i^ the 
magnetic inhomogeneity at the FM/SC interface can in¬ 
duce the non-zero correlations of all three components 
m = 0 and m = ±1 of the spin-triplet state. As a con¬ 
sequence, one appears an extra scattering mechanism, 
called the anomalous Andreev reflection in which elec¬ 
tron incident into the SC is reflected as a hole with the 
same spin (in contrast to the normal Andreev reflection 
in which the incident electron and the reflected hole have 
opposite spins). For the HM/CM/SC junctions, this new 
scattering mechanism, if it exists, leads to the nonzero 
conductance in the subgap region as presented in Fig. [3l 
As one can see the conductance strongly depends on 
the thickness of the conical magnet layer, dcM, i-e. its 
value reaches minimum for dcM = A and maximum for 
dcM = 1.5A, respectively. In Fig. 0] we present the re¬ 
flection and transmission probabilities (i?|J), 

^ function of energy for these two distin¬ 
guished thicknesses. We see that for the CM thickness 
dcM = 1-5A the increase of the conductance in the sub¬ 
gap region is mainly determined by the increase of the 
anomalous Andreev reflection probability On the 

other hand the probability is suppressed for dcM = A 
for which the normal reflection with spin conservation 
(R^l) emerges and also contributes to the conductance 
value leading to its decrease. Regardless of the CM layer 
thickness, for low energy, the anomalous Andreev reflec¬ 
tion probability rJJ^^ drops to zero while the normal re¬ 
flection Rfl increases to unity. This results in the zero 
conductance at zero energy as demonstrated in Fig. [3l 

Now, we discuss the thickness dependence of the con- 
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ductance G{dcM)i important from the viewpoint of ex¬ 
periments in which the critical current is measured as 
a function of the CM layer thickness. In Fig. [5] we 
present the tunneling conductance as a function of en¬ 
ergy and dcM • We see that the conductance oscillates as 



FIG. 5. (Color online) Normalized tunneling conductance G 
as a function of energy E (in units of Ao) and CM layer 
thickness dcM (in units of A). Results for hcM = 1, o = 90" 
and /3o = 0. 

a function of the CM layer thickness with the amplitude 
which is not constant but varies wit dcM- The clear ev¬ 
idence of these oscillations is shown in Fig. [6fb) where 
the cross-section of the G{E,dcM) map is presented for 
the energy E/Aq = 0.02. The small number of points 
in this figure results from the high computational cost of 
the full-self consistent calculations. Results presented in 
Fig.[6](b) are consistent with the phenomena observed in 
experiments, namely (i) the conductance is a nonmono¬ 
tonic function of dcM with peaks for the half-integer 
multiplies of A and (ii) the conductance slowly decays 
with increasing CM layer thickness. For comparison, in 
Fig. [6l^a) the non-self-consistent dependence G{dcM) for 
E/Aq = 0.02 are also presented. As we see the peaks in 
conductance calculated in the non-self-consistent frame¬ 
work are greater than the corresponding peaks calculated 
self-consistently. Moreover, the conductance decay rate 
(with increasing dcM) is slower than in the self-consistent 
approach. 

The irregular nonmonotonic dependence of the conduc¬ 
tance G{dcM) presented in Figs. [5] and [6] can be explained 
as follows. If the spin-active region (the CM layer) is 
present at the FM/SC interface the anomalous Andreev 
reflections can appear giving raise to the nonzero con¬ 
ductance in the subgap regionA^ The probability 
depends mainly in the spin transition in the CM layer. 
Note, that the strength of the spin-flip scattering in the 
CM is proportional to the off-diagonal matrix elements of 
the Hamiltonian ([7]) which have the form —hx{y)Eihz{y). 
Since hx{y) and hz{y) vary periodically, the strength of 
the spin-flip scattering also oscillates with increasing CM 
layer thickness. Nevertheless in the non-adiabatic regime 
the modification of the exchange field seen by electrons 



FIG. 6. (Color online) Normalized tunneling conductance G 
as a function of CM layer thickness dcM (in units of A) for 
energy E/Aq = 0.02. Results obtained by the use of (a) 
the non-self-consistent and (b) the self-consistent approach. 
Results for hcM = 1, o = 90" and /3o = 0. 


flowing through the nanostructure is so fast that the elec¬ 
tronic spins are not able to follow these changes. It en¬ 
tails the irregular oscillations of the anomalous Andreev 
reflection probability and, in consequence, the irregular 
oscillations of the conductance depicted in Figs. [5] and [6l 
In delineating the role of the spin-triplet pairing in the 
charge transport through the HM/CM/SC junctions, it 
is necessary to understand the behavior of the conduc¬ 
tance under the influence of the magnetic configuration 
in the CM layer determined by the value of the exchange 
field amplitude hcM and the angles a and /3o (see Eq.[T]). 
In Fig. [7| we plot the conductance maps G{E,dcM) for 
the exchange field amplitude (a) hcM = 0.2 and (b) 
hcM = 0.3. Other parameters are assumed to be the 
same as in previous calculations, results of which are pre¬ 
sented in Fig. [SI Figure [7] clearly demonstrates the sup- 


(a) hcM=0.2 (b) hcM=0.3 



FIG. 7. (Color online) Normalized tunneling conductance G 
as a function of energy E (in units of Ao) and CM layer thick¬ 
ness dcM (in units of A) for the amplitude of the exchange 
field in the CM layer (a) hcM = 0.2 and (b) hcM = 0.3. 
Other parameters are the same as in previous calculations. 

pression of the conductance in the subgap region with de¬ 
creasing exchange field amplitude hcM • The comparison 
of Fig. Eland Fig.|S]allows to conclude that the strength of 
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this conductance suppression depends on the CM layer 
thickness. It is minimal for dcM = 2.5A. Note that, 
for hcM = 0.2 the conductance peak for dcM = 2.5A 
is still well pronounced [Fig. [7l[a)]. Further reduction of 
the amplitude hcM leads to the situation in which the 
conductance peak survives only for dcM = 2.5A in con¬ 
sistency with the experimental measurements reporting 
the peak of the critical current exactly for this value of 
the CM layer thickness. Although the conductance for 
dcM = 2.5A decreases slower than for other thicknesses, 
even for this value of dcM the conductance in the subgap 
region is suppressed with decreasing hcM- This suppres¬ 
sion is depicted in Fig. [8] which presents G{E) for differ¬ 
ent values of the exchange field amplitude hcM- We see 
that in the limit hcM ^ 0, as expected, the conductance 
in the subgap region tends to zero. We should also no- 



FIG. 8. (Color online) Normalized tunneling conductance G 
as a function of energy E (in units of Ao) for different values of 
the exchange field amplitude hcM- Results for dcM = 2.5A. 

tice clear cups in the conductance for the energy E = Aq. 
The non-self consistent analysis for the FM/SC junctions 
without the CM layer leads to the analytical dependence 
of G(F^ = Ao, hpAd) which is a decreasing function of the 
exchange field in the ferromagnef4^ 

G{E = Ao, hFM) = -— ^ 

As depicted in the insert of Fig. [8] for the HM/CM/SC 
junctions the dependence G{E = Ao, hcM) is an increas¬ 
ing function of hcM in contrast to the FM/SC structure. 
This results from the fact that the conductance for the 
considered energy is mainly determined by the anomalous 
Andreev reflections whose probability increases with in¬ 
creasing amplitude of the spiral magnetic configuration 
in the CM layer. 

The magnetic configuration in the CM layer can be 
modified not only by changing the amplitude hcM but 
also by changing the spatial configuration determined 
by the angles a, /3 and /3o. Figure [9][b) presents the 
conductance G{E) for different angles /3o of a relative 
magnetization of the CM layer measured from the HM 
layer at the HM/CM interface. In panel (a) we present 




FIG. 9. (Golor online) (a) ^-component of the magnetization 
rriz (z) for three different angles f3o of a relative magnetization 
of the GM layer measured from the HM layer at the HM/GM 
interface. Dotted vertical lines marked boundaries of the GM 
layer, (b) Normalized tunneling conductance G as a function 
of energy E (in units of Ao) for different angles /3o. Results 
for dcM = 2.5A, a = 90^ and /3o = 0. 


the 2 ;-component of the magnetization mz{z) for three 
different angles Note that mz{y) is not discontin¬ 
uous but changes smoothly at the interfaces HM/CM 
and CM/SC. It saturates to unity in the HM region on 
the left-hand and penetrates the SC region on the right- 
hand. Although the amplitude of the exchange field h is 
the same for all three cases, the conductance in the sub¬ 
gap region decreases with increasing /3o [see Fig. [9l[b)]. 
This behavior can be easily understood by considering 
two factors. The fist is directly related to the oscilla¬ 
tory dependence of the conductance with the CM layer 
thickness. In fact, the introduction of a relative magne¬ 
tization between the HM and CM layers corresponds to 
the phase shift in the oscillatory dependence of the he¬ 
lical magnetic configuration. Therefore, in the first ap¬ 
proximation, the dependence G{dcM) should be shifted 
in argument by AdcM = 27ra//3o. This shift causes the 
conductance for dcM = 2.5A (corresponds to the max¬ 
imum value for = 0) shifts to lower value related to 
G(2.5A + Ad). The second factor is the increase of the 
normal reflection probability at the HM/CM interface re¬ 
sulting from the discontinuity of the exchange field. The 
presented dependence G{/3o) is important in the simula¬ 
tion of the real HM/CM/SC structure since the conical 
magnets used in the multilayer setup have several ways 
to orient the magnetic moments with respect to the half¬ 
metal magnetization depending on the magnetic coupling 
at the HM/CM interface. 
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All results presented so far have been obtained for a = 
90^ for which the ^-component of the exchange field hy 
is zero - the magnetization in the CM layer rotates in the 
x—z plane. Now, we analyze the conductance for the non¬ 
zero value of hy. In Fig.[T0]we plot the conductance map 
G{E,dcM) for two values of the cone angle (a) a = 30^ 
and (b) a = 60^. We see that for the low cone angle 
a = 30^ the conductance map G{E^dcM) significantly 
differs from the map for a = 90^ (Fig. [5j). Based on 



FIG. 10. (Color online) Normalized tunneling conductance 
G as a function of energy E (in units of Ao) and CM layer 
thickness dcM (in units of A) calculated for (a) a = 30^ and 
(b) a = 60°. Results for hcM = 1, /3o = 0. 


Fig. [To] one can state that regardless of the angle a the 
value of the conductance in subgap region oscillates as a 
function of dcM whereas the period and the amplitude 
of these oscillations are irregular. In Fig. [11] we present 
the conductance map G{E, a) as a function of the cone 
angle calculated for the CM layer thickness (a) d = 2A 
and (b) d = 2.5A nm. As we can see, the conductance in 


B. Adiabatic regime 


In this subsection we analyze the the tunneling con¬ 
ductance in the HM/CM/SC junctions in the adiabatic 
regime for a long period A. In subsection IIII Al we have 
demonstrated that in the non-adiabatic regime the mag¬ 
netization in the CM layer is strongly suppressed com¬ 
pared to the helical exchange field. As presented in 
Fig. [2] (b,c) the amplitude of mx{y) and mz{y) modu¬ 
lation in the CM is about ten times smaller than the 
amplitude of the helical exchange field. Such a strong 
suppression have been obtained for a short period of the 
exchange field modulation A = 3.4 nm corresponding to 
Q <1 (non-adiabatic regime). It has been suggested that 
the magnitude of this suppression can be used as the ad¬ 
ditional parameter to measure the degree of adiabaticity. 



(b) A=5 nm 


rwR 


340 


360 
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380 
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FIG. 11. (Color online) Normalized tunneling conductance G 
as a function of energy E (in units of Ao) and the cone angle 
a for the CM layer thickness (a) d = 2A and (b) d = 2.5A. 
Results for hcM = 1, ydo = 0. 

the subgap region decays with decreasing the cone angle, 
whereas the decay rate strongly depends on dcM- For 
d = 2a it is stronger than for d = 2.5a. 


FIG. 12. (Color online) ^-component of the magnetization 
(red lines) for different periods of the exchange field mod¬ 
ulation (a) A = 3 nm, (b) A = 5 nm, (c) A = 8 nm, (d) 
A = 15 nm. The distributions of the exchange field are plot¬ 
ted by the blue dashed lines, (e) Ratio of the amplitudes 
m'>^axIj^max ^ functiou of thc period of the ex¬ 

change field modulation A. Results for hcM — a — 90°, 
/3o = 0° and dcM — 4A. 

In Fig. [12] we demonstrate the self-consistent z- 
component of the magnetization for different values of 
A assuming hcM = 1- For comparison the distributions 
of the exchange field are plotted by the blue dashed lines. 
As presented in Fig. [12] (compare the red and blue lines), 
the suppression of the magnetization in the CM layer 
is more pronounced for a short period of the exchange 
field modulation A (non-adiabatic regime) and almost 
completely disappears for a long A (adiabatic regime). 
For A = 15 nm there is no difference between the mag¬ 
netization mz{y) and the exchange field hz{y) except 
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FIG. 13. (Color online) Normalized tunneling conductance G as a function of energy E (in units of Ao) and CM layer thickness 
dcM (in units of A) calculated for (a) a = 30^, (b) a = 60^ and (c) a = 90^. Results for A = 15 nm, hcM = 1 and /3o = 0. 


the boundaries of the CM layer where mz{y) smoothly 
changes penetrating the SC region due to the reverse 
proximity effect. Fig. [12] (e) presents the ratio of the 
amplitudes in the CM as a function of the 

exchange field modulation period A. As one can see the 
ratio saturates to the value = i for A grater 

than 15 nm for which the transport through the hetero¬ 
junction can be assumed to be adiabatic. The further 
analysis presented in this subsection will be carried out 
in the adiabatic regime for A = 15 nm corresponding to 
Q^IO. 

Figure [13] shows the normalized conductance as a func¬ 
tion of energy and CM layer thickness for the cone an¬ 
gles a = 30*^, 60^ and 90*^. For a = 90^ corresponding 
to hy = 0 [FigHdTc)] the conductance in the subgap re¬ 
gion oscillates with a period A/2. It reaches maximum for 
dcM being an integer multiple of A/2. In this case (we as¬ 
sume = 0) the spin of electrons injected from the HM 
has the same direction as the exchange field in the CM. In 
the adiabatic regime the spin of electrons flowing through 
the nanostructure follows the exchange field. Therefore, 
in this case the anomalous Andreev reflection probabil¬ 
ity is exactly proportional to the off-diagonal elements 
—hx{y) Eihy{y) which oscillate leading to the regular os¬ 
cillations of the conductance presented in Fig.[T3fc). This 
behavior considerably differs from the irregular conduc¬ 
tance oscillations presented in Fig [5] for non-adiabatic 
regime. Moreover, note, that the value of the conduc¬ 
tance in the subgap region is lower than this obtained in 
the non-adiabatic regime (compare with Fig. [5]) which is 
caused by the enhancement of the normal reflection prob¬ 
ability due to increase of the exchange field modulation 
period A. 

For the cone angle a ^ 90^, corresponding to the non¬ 
zero hy^ the spin of electrons injected from the HM is 
non-collinear with the exchange field at the HM/CM 
interface. Therefore, the electronic spin starts to pro¬ 
cesses around the exchange field direction with the Lar- 
mor frequency which, in the adiabatic regime, is much 
higher than the frequency of the exchange field modula¬ 


tion experienced by electrons flowing through the CM. 
For A = 15 nm (Q ~ 10) one period of the exchange 
field modulation corresponds to ten full-rotations of elec¬ 
tron spin around the exchange field direction. Therefore, 
for the non-zero hy the spin behavior in the CM layer is 
determined by joint effects: Larmor precession and the 
exchange field modulation. For certain energies and the 
CM thicknesses this complex spin behavior leads to the 
enhancement of normal reflection probability with spin 
flip presented in Fig. [Mf b). Note, that the proba¬ 
bility of this scattering mechanism in the non-adiabatic 
regime is close to zero (Fig. [4]). According to Eq. [Tb] in 



23456 23456 

dcJA dew/A 


FIG. 14. (Color online) The reflection and transmission prob¬ 
abilities (a) RJ/, (b) R/J, (c) R/g and (d) T/ as a function 
of energy E (in units of Ao) and CM layer thickness dcM (in 
units of A) for a — 30°. The normal Andreev reflection prob¬ 
ability R/g(E,(icM) = 0. Results for A = 15 nm, hcM — 1 
and ydo = 0. 

the range in which increases, the conductance is sup¬ 
pressed leading to the characteristics G{E,dcM) demon- 
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FIG. 15. (Color online) Normalized tunneling conductance G 
as a function of energy E (in units of Ao) and the cone angle 
a. Results for A = 15 nm and dcM = 4A. 


strated in Fig. m a). In this figure the conductance sup¬ 
pression for the energy above Aq corresponds to green 
areas. Nevertheless this suppression expands also on the 
subgap region for which is even grater than above Aq 
[F ig. [Hl b)]. The evolution of the conductance with in¬ 
creasing cone angle a is presented in Fig. [151 We see that 
the conductance in the subgap region initially decreases 
and then clear conductance peak for a = 30^ appears. 
The position of this peak is shifted on the energy scale 
with increasing cone angle a. 

Finally, in Fig. [16] we present the conductance map 
G(£’, dcM) calculated for a = 90^ and the exchange field 
amplitude dcM =0.2. As shown in previous subsec¬ 
tion in the non-adiabatic regime the decrease of hcM re¬ 
sults in decrease of the conductance in subgap region. 
This conductance suppression is different for different 
CM thicknesses (see Fig. [^ leading to the conductance 
peak for dcM = 2.5A in consistency with the experi¬ 
mental resultsIn the adiabatic regime the exchange 



FIG. 16. (Golor online) Normalized tunneling conductance 
G as a function of energy E (in units of Ao) and GM layer 
thickness (in units of Ao). Results for a — 90" and hcM = 

0 . 2 . 

field amplitude hcM affects the conductance in a dif¬ 
ferent manner. It leads to the conductance decay for 


the thicknesses being an integer multiples of A whereas 
for dcM = {N/2)X the conductance remains almost un¬ 
changed. Therefore, the conductance in the sugap re¬ 
gion oscillates regularly with the period A as presented 
in Fig. [m 


IV. SUMMARY 

We present the detailed analysis of the transport prop¬ 
erties in the HM/CM/SC junction within the fully self- 
consistent framework based on the combined BTK for¬ 
malism and the BdG equations. For comparison, the 
calculations have been also carried out with the use of 
the non-self consistent scheme. One has been shown 
that the peaks in the CM layer thickness dependence of 
the conductance are significantly reduced when the self- 
consistent procedure is applied (cf. Fig. [6]). It is clear 
from our analysis that to properly determine the tunnel¬ 
ing conductance in the HM/CM/SC heterojunctions the 
full self-consistent calculations including the proximity 
effect should be carried out. 

Due to the spin band separation in the HM, the anoma¬ 
lous Andreev reflection mechanism which appears in 
structures with magnetically inhomogeneous layer (such 
as the CM layer) results in the nonzero conductance 
within the subgap region (cf. Fig. [3]). Its probabil¬ 
ity strongly depends on the spin transmission in the CM 
layer. Therefore, we analyze the influence of the exchange 
field modulation A on the behavior of spins of electrons 
transfered through the CM layer. With this respect, we 
show that one can distinguish between two regimes. In 
the non-adiabatic regime (low value of the exchange field 
modulation A) the changes of the exchange field are so 
fast that the spins of electrons are not able to adopt and, 
as a result, they rotate around h irregularly and as an 
average give a small value of magnetization in compari¬ 
son to the amplitude of the exchange field. On the other 
hand, in the adiabatic regime (high value of A) the mag¬ 
netization coming from the spins of electrons is almost 
identical to the exchange field position dependence within 
the CM layer (cf. Fig [12]). The conductance behavior as 
function of the energy and the CM layer thickness has 
different behavior in the two mentioned regimes (cf Figs. 
[5land [T3lf c)). The regular oscillations observed in the adi¬ 
abatic regime are caused by the fact that anomalous An¬ 
dreev reflection probability is exactly proportional to the 
off-diagonal elements of the Hamiltonian —hx{y)-\-ihy{y). 
This leads to constant hight of the peaks in the CM layer 
thickness conductance dependence, whereas for the case 
of the non-adiabatic regime the hight of those peaks is 
decreasing with increasing thickness, dcM^ in consistency 
with the experiment observation (cf. Fig. [6]). Moreover 
in the non-adiabatic regime the decrease of the exchange 
filed amplitude results in the well pronounced conduc¬ 
tance peak for dcM = 2.5A for which the peak of the 
critical current was observed in the experiment The in¬ 
fluence of other parameters characterizing the exchange 
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field behavior in the CM layer are also analyzed (such as 
the a, P and /3o angles). It is shown that the conduc¬ 


tance is strongly affected by the geometrical structure of 
the exchange field determined by the cone angle a and 
the rotational a4.21a (PWD, AO, DPC) hacked 
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